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ABSTRACT 

Context. The cusp-core discrepancy is one of the major problems in astrophysics. It results from comparing the observed mass 

distribution of galaxies with the predictions of Cold Dark Matter simulations. The latter predict a cuspy density profile in the inner 

parts of galaxies, whereas observations of dwarf and low surface brightness galaxies show a constant density core. 

Aims. We want to determine the shape of the dark matter potential in the nuclear regions of a sample of six nearby irregular dwarf 

galaxies. 

Methods. In order to quantify the amount of non-circular motions which could potentially affect a mass decomposition, we first 
perform a harmonic decomposition of the Hi Hermite velocity fields of all sample galaxies. We then decompose the Hi rotation 
curves into different mass components by fitting NFW and pseudo-isothermal halo models to the H i rotation curves using a 
minimisation. We model the minimum-disc, the minimum-disc + gas, and the maximum-disc cases. 

Results. The non-circular motions are in all cases studied here of the order of only a few km s"' (generally corresponding to less than 
25% of the local rotation velocity), which means that they do not significantly affect the rotation curves. The observed rotation curves 
can better be described by the cored pseudo-isothermal halo than by the NFW halo. The slopes of the dark matter density profiles 
confirm this result and are in good agreement with previous studies. The quality of the fits can often be improved when including the 
baryons, which suggests that they contribute significantly to the inner part of the density profile of dwarf galaxies. 

Key words, dark matter - galaxies: dwarf - galaxies: kinematics and dynamics 



1. Introduction 

Early observations of spiral galaxies have shown that their dy- 
namical mass, as inferred from their surprisingly flat rotation 
curves, cannot be explained by luminous matter alone (e.g., 
|Bosma| 1978) [Rubin et al.|1978| l. Since then, "dark matter" has 
become part of modern cosmology. Currently, the most suc- 
cessful models use a cosmological constant A as well as colli- 
sionless and dissipationless Cold Dark Matter (CDM). ACDM 
simulations like, e.g., the Millennium Simulation by Springel 



|et al. ( |2005| l have been very successful in describing the ob- 



served large-scale structures in the Universe ( Spergel et al. 2003 
|Spergel et al.||2007) l, but they do not work properly on galaxy 
scales. 



Next to the "missing satellite" problem ([Moore et al.|1999) , 
the "cusp-core" discrepancy (e.g., [de Blok et al.[[2001b| l is still 
causing a lot of debate between observers and cosmologists. 
Numerical simulations predict cuspy haloes in the inner parts 
with a density distribution described by a power law p{r) ~ r" 
with a ranging from -1 (e.g., [ Navarro et al.][1996 1 to -1.5 (e.g., 
[Moore et al.| 1998, 1999). This cusp leads to a steeply rising 
rotation curve. However, observations of dwarf and low sur- 
face brightness galaxies show that their rotation curves rise less 
steeply than predicted by CDM simulations (e.g., |de Blok"&| 
[Bosma[[2002[ l. Similar conclusions have also been made from 
observations of high surface brightness disc galaxies (Salucci 
[200 l| l. At small radii (typically few kpcs), the mass distribu- 



tion can better be described by a central, constant-density core 
(e.g., Fl ores & Primack|1994{[C6te et al. 2000; Marchesin i et al. 



This interpretation has been met with skepticism by the cos- 
mologists. Thus, systematic eff'ects in the data like beam smear- 
ing, slit misplacement, slit width, and seeing, as well as the 
inclination of the galaxy and non-circular motions have been 
used to argue against the fact that the observations are incom- 
patibl e with cusps ([van den Bo sch et al.||2000[ [Swaters et al. 



2004; Spek kens et al.[[2005| l. [Rhee et al 



2003| [Hayashi et al. 

(2004) tried to quantify several of these errors and found that 
they add up to 30% - 50% underestimation of the rotation curve. 
However, repeated one-dimensional long-slit spectra observed 
by independent observers at different telescopes ( de Blok et al. 
2003 ) as well as high-resolution Ha observations using a two- 
dimensional velocity field ( Simon et al.[j2005 Kuzio de Naray 



et al.i2006„2008„Spano et a 



2008 1 rule out uncertainties related 



to slit spectroscopy and still suggest the presence of an isother- 
mal dark matter core in the inner parts of disc galaxies. Non- 
circular motions can indeed affect the results of rotation curve 
studies as it is generally assumed that the particles are on cir- 
cular orbits. Several authors report that CDM haloes are triaxial 



objects with a globally elongated potential (e.g., Hayashi et al. 
[2004 ' Hayash i & NavarropO Oe*). They suggest that this triaxial- 
ity induces large non-circular motions in the inner parts of galax- 
ies (up to 15% of the maximum rotation velocity), leading to the 
observed cored profiles. However, [de Blok et al.[ ( |2003| l simu- 
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lated rotation curves showing that non-circular motions of the 
order of 20 km s ' over a large fraction of the disc are needed to 
make them consistent with CDM haloes. Furthermore, recent ob- 
servations by [Gentile et al. (2005 ) and Trachter nach et al.| ( |2008| l 
reveal that non-circular motions are typically of the order of few 
kms"', i.e., too small to explain the cusp-core discrepancy. This 
might not always be the case as |Spekkens & S ellwood' ( |2007| l 
found high non-circular motions in the spiral galaxy NGC 2976, 
which have a major effect on the rotation curve and which they 
explained by a bar in the inner 500 pc. 

In this paper, the cusp-core discrepancy is addressed by us- 
ing the rotation curves derived from H i synthesis data for a mass 
decomposition of a sample of six nearby irregular dwarf galax- 
ies. The galaxies have been chosen to be a subsample of dwarf 
irregular galaxies in the Local Volume (D < lOMpc), which 
have been observed in H i with sufficiently high spatial resolu- 
tion (< 1 kpc). Some general properties are given in Table[T] 

In order to rule out systematic effects in the data, we first 
measure the non-circular motions. Non-circular motions can 
have two major causes: chaotic non-circular motions can be in- 
duced for example by star formation (Oh et al. 2008), systematic 
non-circular motions relate to the potential (e.g., spiral arms, tri- 
axiality of the halo, see Schoenmakers et al.|1997 1. By perform- 
ing a harmonic decomposition of the velocity fields, we are able 
to quantify the systematic non-circular motions. 

This paper is organised as follows: in Sect. 2, the different 
data sets and the creation of the velocity fields are briefly de- 
scribed. The velocity fields and rotation curves are presented in 
Sect. 3. In Sect. 4, we describe the procedure of the harmonic 
decomposition and subsequently show and discuss our results. 
Section 5 presents the theoretical background of the mass mod- 
elling, which is followed by a discussion of the results from the 
mass decomposition in Sect. 6. In Sect. 7, the main results are 
summarised. 



2. Observations and data reduction 

For this study we selected six nearby irregular dwarf galaxies. 
Some observational details are given in Table [2] 

Hi data for the two northern galaxies, NGC 2366 and 
NGC 4861, were obtained with the Very Large AiTay (VLA). 
The IB(s)m galaxy NGC 2366 forms part of "The Hi Nearby 

Its peculiar 



Hermite h3 polynomials to all line profiles using the GIPSY | 



(The Groningen Image Processing System; van der Hulst et al. 
1992[ ) task xgaufit. This method allows to accurately define the 



Galaxy Survey" (THINGS, [Walter et al.][2008 



kinematics were the subject of several recent publications Jde, 
Blok et al.|[200 8 ; Trac hternach eraD 2008 , Oh ^eFarjlSOOSl [vaij 
Eymeren etal.|2009 a, hereafter dB08, CT08, Oh0 8, and vE09a). 
The S B(s) m galaxy NGC 486 1 was studied by [Wilcots et al. 
1996| l and Thuan et al. ( 2004|i. A detailed kinematic analysis of 



this galaxy by van Eymeren et al. ( 2009b i is based on the com- 
bined H I data of two C array and one D array observation. The 
data reduction process of NGC 2366 is described in [Walter et al.| 
( [2008.) , the one of NGC 4861 in van Eymeren et al.[ ( [2009b l 

The four southern dwaif galaxies, ESO059-G001, ESQ 215- 
G?009, NGC 5408, and IC5152, are selected from flie "Local 
Volume Hi Survey" (LVHIS0 [Koribalski & Jerjen 



2008 



Koribalski et al.|2009[ l. H i data for these galaxies were obtained 



with the Australia Telescope Compact Array (ATCA) and re- 
duced by the LVHIS team (for details see [Koribalski et al.[2009l l. 
All galaxies were observed for 12 hours each in at least three dif- 
ferent configurations (see Table |2]). 

For the purpose of this paper Hi velocity fields are cre- 
ated from the naturally-weighted data cubes by fitting Gauss- 



peak velocities, whereas the standard intensity-weighted mean 
velocities are biased toward the longest tail of the velocity pro- 
files as soon as the distribution is not symmetric (see also dis- 
cussion in dB08). In order to separate true emission from noise, 
we set some limits for the fitting algorithm (see dB08): the fit- 
ted profiles need to have amplitudes higher than icrchan, where 
Cchan IS the average noise in the line-free velocity channels of the 
cube. The minimum dispersion to be fitted has to be higher than 
the spectral resolution of the cube (see Table |2]for the noise val- 
ues and channel separations). After creating the velocity maps, 
a small amount of noise pixels has to be removed that was ad- 
mitted by the filter criteria. Therefore, we follow again dB08 and 
use the integrated H i column density map as a mask: fits are only 
retained if the total flux in the integrated H i map is higher than 
3a- where ctn is the noise in the integrated H i map. It is defined 
as 's/NcTchaiu where is the number of channels with signal con- 
tributing to each pixel. No smoothing is applied. Altogether, this 
gives us high quality velocity fields. 

3. Tlie velocity fields and the rotation curves 

Figure[T]shows the H i velocity fields (left panels). In order to de- 
rive a rotation curve, we use the GIPSY task rotcur: tilted rings 
of half the beam width are fitted to both sides of the H i veloc- 
ity fields. The systemic velocity, the centre position, the inclina- 
tion, and the position angle of each ring are iteratively defined 
by keeping all parameters fixed except for the one we want to 
measure. As initial estimates we use kinematic parameters de- 
rived from fitting isophotes to the integrated H i intensity map 
(GIPSY task ellfit). The final rotation curve is then created by 
keeping the best-fitting parameters fixed. We also derive rota- 
tion curves for receding and approaching side only by keeping 
the systemic velocity and the centre position fixed to the values 
obtained in the joint approach and by defining the inclination 
and the position angle iteratively. 

The resulting rotation curves of all sample galaxies are 
shown in the middle panels of Fig. [T] The error bars of the ro- 
tation curves indicate the rotation curves of receding (top) and 
approaching side (bottom). We then take the iteratively defined 
kinematic parameters from the tilted-ring analysis (see Table [3]l 
in order to create a model velocity field for each galaxy, which 
is subsequently subtracted from the observed velocity field. This 
shows how well the kinematics of each galaxy have been mea- 
sured in the tilted-ring analysis. The residual maps are presented 
in Fig. [T] right panels. In most of the cases, the derived pa- 
rameters describe the kinematics of the galaxies quite well and 
the absolute values of the residuals are below an absolute value 
of lOkms A comparison of the positions of the optical and 
dynamic centres shows that the offsets (given in the last two 
columns of Table |3]l are generally smaller than one beam size. 

A different test to prove the quality of the rotation curves 
was done, e.g., in vE09a for NGC 2366. They chose different 
approaches in order to derive the rotation curve. First, the initial 
estimates from ellfit were kept fixed. Then, a rotation curve with 
the best-fitting parameters kept fixed was derived (see above). As 
a third approach, the best-fitting parameters were left free. vE09a 
could show that the deviations between these three approaches 
are generally small with higher deviations in the outer parts of 
the rotation curve. In the inner 1 kpc, the deviations are less than 



LVHIS project webpage: www.atnf.csiro.au/research/LVHIS/ 
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Galaxy 


Type 


optical centre 


D 






a (J2000.0) 


6 (J2000.0) 


[Mpc] 


(1) 


(2) 


(2) 


(2) 


(3) 


NGC 2366 


IB(s)m 


O?*" 28™ 54.6^ 


+69° 12' 57" 


3.44 


ESO059-G001 


IB(s)m 


O?*" 31™ 18.2-^ 


-68° 11' 17" 


4.57 


ESO215-G7009 


SAB(s)m 


IQh 57m 29.9^ 


^8° 10' 43" 


5.25 


NGC 4861 


SB(s)m 


12h 59™ 02.3-^ 


+34° 51' 34" 


7.50 


NGC 5408 


IB(s)m 


14'' 03™ 20.9^ 


^1° 22' 40" 


4.81 


IC5152 


IA(s)m 


22h 02™ 41.5-^ 


-51° 17' 47" 


2.07 



Notes : (1) th e name of the galaxy; (2) data from NED; (3) distance ref erences: NGC 2366: [Tolstoy et al.] ( (T995l l, ESQ 059-GOOl : |Karachentsev| 
leTaTlpOOg), ESO215-G7 009: K arachentsev et al.| pOPTl l, NGC4861: |de Vaucouleurs et al.| ( |1991| irNGC5408: [Karachentsev et al.M2002a| i7 
IC 5152: |Karachentsev et al.H2002bi r 

Table 2. The main observational parameters. 



Galaxy 


Telescope 


Arrays 


Spectral resolution 


Synthesised beam 


Noise 


References 








[kms-'] 


" X " 


pc x pc 


mjy beam ' 




NGC 2366 


VLA 


B+C+D 


2.6 


13x 12 


217 X 200 


0.5 


(1) 


ESO059-G001 


ATCA 


EW352+750C+1.5B 


4 


51x47 


1130X 1041 


1.3 


(2) 


ESO215-G7009 


ATCA 


EW352+750A+1.5C+6A 


4 


21x 21 


530x 530 


1.0 


(2) 


NGC 4861 


VLA 


C+D 


5.2 


31x 30 


1127X 1091 


0.6 


(3) 


NGC 5408 


ATCA 


375+750D+1.5A 


4 


57x 50 


1329 X 1166 


1.9 


(4) 


IC5152 


ATCA 


EW367+750A+1.5A 


4 


50x43 


502 X 432 


1.5 


(2) 



References: (l)|Walter et aTlpOOSll; (2)|Koribalski et al.|p009ll; (3)|van Eymeren et al.|([2009b); (4)|van EymerenlpOOS). 



+2kms We perform this test for all galaxies of our sample 
and generally measure equally small deviations in velocity, at 
least in the inner kpc. 

The neutral gas distribution in ESO215-G7009 and its kine- 
matics (including the derivation of a rotation curve) have been 
studied in detail by Warren et al. (2004). Rotation curves of 
ESQ 059-GOOl and ESO215-G7009 have also been derived by 
|Kirby et"ar| ( |2009] ). A comparison reveals that their parameters 
are in good agreement with our values. 



Subsequently, a model velocity field is created and subtracted 
from the original velocity field producing a residual velocity 
field. Note that the centre coordinates should be kept fixed (see 
S ch oenmakers 1 1 999\ . We run the routine with two different pa- 
rameter sets: in the first case, all parameters except for the centre 
position (listed in Table[3]) are left free. In the second case, incli- 
nation and position angle are fixed to the values given in Table[3] 
The quadratically-added amplitude for each order of the har- 
monic decomposition is calculated using 



4. Harmonic decomposition 

In order to search for non-circular motions, we decompose the 
velocities detected along the tilted rings into multiple terms of 



sine and cosine. Following Schoenmakers ( 1999 1, the line of 
sight velocity vjos can be described as: 



Vlo: 



;('") = v'sys(r) + ^ c,n(r) COS m ijj + s^{r) sin m 



(1) 



m— 1 



where is the maximum fit order used, r is the radial distance 
from the dynamic centre, i/r is the azimuthal angle in the plane 
of the disc, and Vjys is the 0th order harmonic component cq. 
As shown in CT08, a decomposition of the velocity field up to 
third order is sufficient to capture most of the non-circular signal. 
Therefore, we restrict our analysis to = 3. In case of purely 
circular motion, only m - and m - 1 terms are included into 
Eq.g 

The decomposition is performed with the GIPSY task reswri, 
which fits a tilted-ring model to the original velocity field us- 
ing circular rotation, and decomposes the line-of-sight veloc- 
ity along each ring into multiple terms of sine and cosine. 



Ai(r) = (2) 
for m - 1 (ci is the circular rotation velocity) and 



A™(r) 



ctnir) + 4(r) 



(3) 



for m > 1. Additionally, the quadratically-added amplitude of 
all non-circular harmonic components (up to = 3 in this case) 
was calculated from 



AM 



4 



s-ir) + c^ir) + i^(r) + c^(r) + i^(r). 



(4) 



We also derive the median value of A,„(r) for each harmonic or- 
der m in two different ways: once, using the entire radial range 
and once, using the inner 1 kpc only, which is the region where 
the distinction between a cusp and a core becomes most obvious 
(de Blok 2004|). More information about harmonic decomposi- 
tions in general can be found in Schoenmakers et al. ( 1997 j ). 

In the following subsections the results for each galaxy are 
discussed separately. For a better understanding of the interpre- 
tation see the rules of thumb given in |Schoenmakers| ( |l999| l. 
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22^ 3m 5 10 22^ 3"" 2 

R.A. (2000.0) r (kpc) R.A. (2000.0) 



Fig. 1. The H i kinematics of all sample galaxies. Left panels: Hermite h3 velocity fields. The systemic velocity is marked by a bold 
line, the beam is placed into the lower left comer Contour levels are overlaid in steps of 10km s '. Middle panels: the rotation 
curves derived from a tilted-ring analysis (see Table[3]for the fit parameters). Receding and approaching side are indicated by error 
bai-s. Right panels: the residual maps after subtracting a model velocity field (created from the fit parameters in Table [3]l from the 
observed velocity field. 
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Table 3. Kinematic parameters of the sample galaxies. 
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Galaxy 


dynamic centre 




V,ot 


i 


PA 


'^opt ^dyn 






a (J2000.0) 


6 (J2000.0) 


n 1 -1 

[kms '] 


[kms '] 


n 


[°] 


["] 


["] 


(1) 


(2) 


(2) 


(2) 


(2) 


(2) 


(2) 


(3) 


(3) 


NGC 2366 


07'" 28™ 53.6'' 


4-69° 12' 28" 


98 


50 


63 


43 


-5 


-29 


ESO059-G001 


07'' 31"' 18.0^* 


-68° 11' 13" 


528 


69 


44 


325 


-1 


+2 


ESO215-G7009 


10'' 57"' 33.0' 


^8° 10' 47" 


599 


11 


28 


119 


+31 


-4 


NGC 4861 


12'' 59™ 01.4^* 


+34° 51' 44" 


835 


46 


67 


15 


-11 


+ 10 


NGC 5408 


14'' 03"' 22.7^* 


^1° 22' 12" 


502 


53 


58 


300 


+20 


+28 


IC5152 


22'' 02™ 41.9'' 


-51° 17' 45" 


123 


44 


52 


284 


+4 


+2 



Notes: (1) the name of the galaxy; (2) kinematic parameters derived by fitting a tilted-ring model to the Hi Hermite velocity fields, the position 
angle is measured counter-clockwise from north to the receding side of the galaxies; (3) position offset of the dynamical centre in comparison to 
the optical centre. 



4.1. NGC 2366 

Using accurately defined centre coordinates, CT08 have aheady 
performed a harmonic decomposition of the Hermite velocity 
field of NGC 2366 (see also dB08). We also use the Hermite ve- 
locity field, but different centre coordinates (see Table [3]l. The 
radial distribution of all fitted parameters is shown in FigT|2] The 
black filled triangles represent the case where all parameters ex- 
cept for the centre position are left free. The inclination varies 
significantly in the inner 1.5 kpc. As NGC 2366 is dominated 
by close-to solid body rotation, it is difficult to determine Vrot 
and i simultaneously. At a radius of 1 kpc, cq rises from 95 to 
102 kms"' and the position angle jumps from about 35° to 55°. 
At the same radius, the inclination has very low values, which 
leads to very high values of the inclination corrected c\. This 
is probably caused by the sudden decrease in H i intensity. At a 
radius of about 3 kpc, the si and 53 terms show the character- 
istic wiggles caused by spiral arms. As discussed by Tikhonov' 
[&Galazutdinova (2008 ) and vE09a, NGC 2366 shows evidence 
for two weak spiral arms. 

The light grey crosses in Fig. |2] show the constrained case 
where we use a fixed inclination of 63° and a fixed position angle 
of 43° (see Table |3]l. For most of the tilted rings the values are 
in very good agreement with the ones of the unconstrained fit. 
A pronounced deviation can be seen in the s\ term where the 
constrained fit leads to higher amplitudes between 1 and 2 kpc 
than the unconstrained fit. This is due to the constant position 
angle which differs significantly from the unconstrained values. 

The median value of the quadratically-added amplitude for 
each order is given in the lower left panel of Fig. [2] In both cases 
the amplitudes are below 3kms As the innermost region of 
a galaxy is most significant for the cusp-core debate (de Blok 
|2004) , we also show the median values within 1 kpc (open tri- 
angles). These are generally below 1 kms ' and therefore even 
smaller than the ones averaged over the entire radial range. A^ir) 
is below 3 kms"' in the inner 2.5 kpc and rises to 9 kms"' in 
the outer parts contributing about 9% (constrained case) and 
6% (unconstrained case) to the maximum rotation velocity of 
NGC 2366 (see Table|4]l. 

In general, our results agree very well with the results of 
CT08 despite the different centre coordinates. This is to be ex- 
pected as CT08 have already addressed the problem of inaccu- 
rate centre positions. They showed that for dwarf galaxies, which 
are characterised by solid body rotation in the inner parts, the 
results of the harmonic decomposition are insensitive to small 
offsets in the centre position. 



4.2. ESO059-G001 

The results of the harmonic decomposition of ESO059-G001 
are shown in Fig. [3] The inclination at 2 kpc drops to 30°, which 
leads to deviations of the rotation velocity. The position angle 
varies only slightly over the entire radial range. The amplitudes 
of the non-circular components are the smallest of all galax- 
ies studied here. |Parodi et al.| ( |2002| l report "rudiments of spiral 
arms" in ESO 059-GOOl, which can, however, not be seen in the 
harmonic components (probably due to their weakness and the 
relatively low spatial resolution of the Hi data). 

The median amplitudes are below 1 km s"' . Aj(r) stays below 
3 kms"' over the entire radial range and only contributes 2% 
(constrained case) and 1 % (unconstrained case) to the maximum 
rotation velocity (see Table |4|i. 

ESO 059-GOOl is dominated by close-to solid body rotation 
and additionally, the inclination is low. A tilted-ring analysis 
with free parameters is therefore difficult to perform. However, 
the constrained case with a fixed inclination of 44° and a posi- 
tion angle of 325° still gives acceptable values. As long as / does 
not differ too much from the median value, the fits are in good 
agreement with the unconstrained ones. The residual map given 
in Fig. [T] also shows that the kinematic parameters derived from 
the constrained approach represent the galaxy quite well. 

4.3. ESO215-G7009 

ESO215-G7009 has a strongly varying and very low inclination 
in the inner 3.5 kpc (see Fig. [4]|. As shown for ESO 059-GOOl, 
a low inclination makes it very difficult to perform a tilted-ring 
analysis by leaving all parameters free. The partially unreason- 
able values for the inclination in the inner 3.5 kpc of the un- 
constrained case (with / as low as 5°) lead to an unphysical 
rotation curve. A constrained fit with fixed inclination and po- 
sition angle (see Table [3]) gives much more reasonable results 
which are, apart from the inner 3.5 kpc, largely consistent with 
the results from the unconstrained case. Therefore, we decide 
to ignore the unconstrained case for the subsequent analysis of 
ESO215-G7009 and to use the constrained case instead. 

Most of the amplitudes of the harmonic components are 
close to zero. The values of Ay{r) are typically below 10 kms"' 
contributing about 8.5% to the maximum rotation velocity (see 
Table[4]i. The median amplitudes are below 5 km s"' , in the inner 
kpc below 3 km s"' . 

We also show the results of the unconstrained fit: the low 
inclination causes very large circular velocity components (note 
that we do not plot the values between 1 and 3 kpc that have 
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Fig. 2. The results of the unconstrained (black triangles) and constrained (light grey crosses) harmonic decomposition of NGC 2366. 
For the constrained fit the inclination and the position angle were fixed to the values given in Table [5] Upper left: following Eq.[T| 
circular (ci) and non-circular (c2, cj, si, S2, s^) harmonic components corrected for inclination and plotted vs. radius r. Upper 
right: systemic velocity cq, inclination /, and position angle PA, again plotted vs. radius r. The dashed horizontal lines represent 
the error weighted means of the unconstrained fit. The inclination and position angle are from the tilted-ring fit assuming circular 
rotation. For the unconstrained case the error bars represent the formal errors. Lower left: the median amplitudes of the individual 
harmonic components derived from Eqs.|2]and[3]plotted vs. harmonic number m. The error bars denote the upper and lower quartile 
of the distribution of the unconstrained < >■ The open triangles represent the median amplitudes within the inner 1 kpc. Lower 
middle: Ar(r), the quadratically-added amplitude of all non-circular components derived from Eq.|4] The uncertainties are estimated 
using Gaussian error propagation. Lower right: the H i intensity distribution. The dynamic centre is marked by a white cross. 



high formal errors). The low inclination becomes also notice- 4.4. NGC 4861 
able in the non-circular components, which are partly as high as 

20kms The results of the harmonic decomposition of NGC 4861 are 

given in Fig. |5] The amplitudes of the non-circular components 
are usually close to zero. No pronounced deviation can be seen, 
which fits into the picture of a quiescent galaxy without spiral 
arm structure as described in van Eymeren et al. ( 2009b i. The 
median values of the quadratically-added amplitudes for each 
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Fig. 3. The same as in Fig.|2]for ESO059-G001. The inclination and the position angle for the constrained fit are listed in Table [3] 



harmonic order are below 2.5 km s"' over the entire radial range 
and even below 1 kms ' in the inner 1 kpc. The distribution of 
A,(r) shows that the amplitude is generally below 3 kms"', ex- 
cept for the outermost part. On average, the non-circular motions 
contribute about 6% (constrained case) and 3% (unconstrained 
case) to the maximum rotation velocity. The amplitudes of the 
harmonic components of the constrained fit (inclination of 67°, 
position angle of 15°) usually agree to within 2 kms"' with the 
ones from the unconstrained fit. 



fit. Nevertheless, the amplitudes of the non-circular motions are 
small. Ar(r) is below 3 kms"' in the inner 5 kpc and rises to 
about 6 kms"' in the outer parts. On average, the non-circular 
motions contribute less than 8% (constrained case) and 4% (un- 
constrained case) to the maximum rotation velocity. The median 
amplitudes of the individual harmonic orders are below 3 km s"' 
independent of the radial range we look at. 



4.5. NGC5408 

Figure|6lshows the results for NGC 5408. In her PhD thesis, 'van' 
[Eymeren (2008) suggested that this galaxy has two kinematic 
systems, which is probably the reason for the clear and well- 
defined radial variation of the inclination in the unconstrained 



For the constrained case we use a fixed inclination of 58° and 
a fixed position angle of 300° (see Table[3]l. In the inner 2.5 kpc 
this inclination agrees with the one from the unconstrained fit. 
However, it differs significantly from the inclination of the outer 
kinematic system. For r > 2.5 kpc this results in different in- 
clination corrections for the constrained and the unconstrained 
cases, which can be seen, e.g., in the small offset in ci. 
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Fig. 4. The same as in Fig.|2]for ESO 215-G7009. The incUnation and the position angle for the constrained fit are Usted in Table [3] 
The analysis of this galaxy is based on the constrained case (see Sect. |4.3[ ). Therefore, no < > and Ai(r) are given for the 
unconstrained case. 



4.6. IC5152 

Except for one outlier in the inner 1 kpc, the inclination of 
ICS 152 is roughly constant with radius (see Fig. Ij. The posi- 
tion angle increases at 2 kpc by about 20°. Altogether, this leads 
to small non-circular components, often close to 0km s '. The 
constrained case (with an inclination of 45° and a position an- 
gle of 287°) and the unconstrained case are in good agreement, 
although the fixed position angle leads to higher amplitudes in 
the ii term. Ar(r) varies over the entire radial range, but stays 
below 8km s '. The difference between the amplitudes of the 
constrained and unconstrained values are several kms which 
is due to the difference in the si term. The median amplitudes 
of the individual harmonic orders are below 2 km s ' for the un- 
constrained case. In order to calculate the median amplitudes in 



the inner 1 kpc, we removed the outlier at 0.5 kpc. Especially the 
m - \ term of the constrained case is high, which is again due 
to the high amplitudes of the s\ term. In comparison to the other 
galaxies of our sample, IC5152 shows the highest contribution 
of non-circular motions to the maximum rotation velocity. 

4.7. Discussion 

In all sample galaxies, the non-circular motions contribute on av- 
erage about 8% (constrained case) and about 4% (unconstrained 
case) to the maximum rotation velocity. However, in the in- 
ner few kpcs the rotation velocity has not reached its maxi- 
mum value, yet. Furthermore, for the cusp-core debate the in- 
ner few kpcs are most relevant, whereas non-circular motions 
in the outer parts of a galaxy, though surely interesting, are of 
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Fig. 5. The same as in Fig.|2]for NGC4861. The inclination and the position angle for the constrained fit are listed in Table [3] 



no concern. Therefore, we perform a more detailed analysis of 
the quadratically-added amplitude of all non-circular motions 
within the inner 5 kpc. We concentrate on the results of the con- 
strained case. Especially in the inner few kpcs, a strong change 
in inclination or position angle is not to be expected so that it is a 
reasonable assumption to keep both parameters fixed. As Table[4] 
shows, the amplitudes of the constrained case are usually up to 
a factor of two higher than the amplitudes of the unconstrained 
case. Therefore, the constrained case gives some kind of upper 
limit to the amount of non-circular motions. 

Figure |8] shows the amplitudes of the constrained case aver- 
aged within rings of 1 kpc width (left panels) and within rings 
of increasing radius (right panels). Furthermore, we differen- 
tiate between the absolute amplitudes (upper panels) and the 
amplitudes which are normalised by the local rotation velocity 
(lower panels). In all sample galaxies and for all radii, the non- 
circular motions generally contribute less than 25 % to the local 



rotation velocity, sometimes even less than 10%. In the inner 
1 kpc, < A,(r) > is often close to 1 or 2kms"' (see Table |4]l, 
which corresponds to about 5 to 25% of the local rotation veloc- 
ity. According to the simulations by jHayashi et al.| ( |2004 ) and 
IHayashi & Navarro ( |2006| l, the non-circular motions add up to 
about 50% of the local rotation velocity at a radius of 1 kpc and 
even more below 1 kpc. 

In order to check if our harmonic decomposition is able to 
quantify most of the non-circular motions, we follow CT08 and 
make use of the residual velocity fields. In Fig.|9]we compare 
the median of the absolute residual velocity fields created with 
reswri (constrained case) with the median of the absolute resid- 
ual velocity fields created with rotcur. The galaxies from our 
sample are presented by open black triangles. For a comparison 
we also show the data points from CT08 (filled light grey tri- 
angles). In agreement with CT08, the rotcur residuals are larger 
than the reswri ones as rotcur does not take non-circular motions 
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Fig. 6. The same as in Fig.|2]for NGC 5408. The inclination and the position angle for the constrained fit are listed in Table [3] 



into account. In comparison to most of the galaxies studied by 
CT08, our sample dwarf galaxies lie at the lower end of the dis- 
tribution. The median value of the reswri residuals < Mres > is 
below 2kms ' (see Table |4|i, which implies that the harmonic 
decomposition up to third order has captured most of the non- 
circular motions. 



Altogether, we could show that while large non-circular mo- 
tions might be found in dwarf galaxies, they are the exception 
rather than the rule. The galaxies in our sample have non-circular 
motions that are that small that they do not significantly affect 
the rotation curves. This implies that non-circular motions can- 
not artificially flatten the slope of the density profile and turn a 
cusp into a core. 



5. Mass models 

The Hi rotation curves of our sample galaxies are only 
marginally influenced by non-circular motions. Therefore, we 
use the ones derived with rotcur (see Sect. |3]and Fig.[r| middle 
panels) to perform a mass decomposition with the GIPSY task 
rotmas. This task allows to interactively fit the different compo- 
nents (halo, gas, stars) to the observed rotation curve by min- 
imising the of the parameter space. In this section, the dark 
matter halo models are introduced and the derivation of the mass 
components is described. 

5.1. Dark matter halo models 

We use two models to describe the dark matter halo, which are 
the cuspy NFW halo (e.g., [Navarro et al.||1996 i and the cored 
pseudo-isothermal (ISO) halo (e.g., |Binney & Tremaine|1987| l. 
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Numerical simulations show that the density of CDM haloes 
rises steeply towards the halo centre (e.g., Navarro et al., 1996) 1. 
The NFW mass-density distribution is described as 

Pi 



^2' 



(5) 



r/rs ■ (1 + r/r,r 

where p; is related to the density of the Universe at the time of 
halo collapse, and is the characteristic radius of the halo. This 



leads to a rotation curve of 

Vrot('") = V'200 



ln(l + cx) - (cx)/( l + cx) 



^ x[ln(l -Hc)-c/(l H-c)] 



(6) 



with X = r/r2oo- The concentration parameter c - r2oo/'"s is di- 
rectly correlated to V200, the circular velocity at r2oo <de Blok 
let al.|2003| l 



log c = 1 . 1 9 1 - 0.064 log V20() - 0.032 log v: 



200' 



(7) 



where r2oo is the radius at which the density contrast exceeds 
200, i.e., roughly the virial radius (Nava rro et al.| [T996). Note 
that we do not correct the NFW fits for adiabatic compression. 
As|Sellwood & McGaugh (2005 1 show, the central halo concen- 
tration of dwarf galaxies increases only slightly from the pri- 
mordial to the compressed halo, which means that it still fits the 
result from N-body simulations (a = 1). 

As a second model the empirically derived ISO halo is used 
( Binney & Tremaine|1987| l. It describes a dark matter halo that 
has a core of roughly constant density. The density profile is 
given by 



PisoC'") = po 



(8) 
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Table 4. Derived quantities from the harmonic decomposition. 



Galaxy 
(1) 


< A, > 

[Km S J 

(2) 


< A,..ikpc > 
[Km S J 
(3) 


< A, >/v,.o, 
L%J 
(4) 


< M,.es > 
[Km S J 

(5) 


< A,,, > 

[Km S J 

(6) 


< A,,,.ikpc > 

[%\ 

(7) 


< Ar >A'iot 
[Km S J 

(8) 


< Mes.c > 

LkpcJ 
(9) 


'"max 

(10) 


NGC 2366 
ESO059-G001 
ESO215-G7009 
NGC 4861 


9 Q 1+4.09 
^■^^-1.06 

' -0.09 


0.77;:: 

0.78;:: 


5.8 
1.1 

3.3 


2.15 
0.43 

1.04 


4 81+2 ^2 

1 4^+0.62 

6 i7+i M 

U.J /_j 38 

2 60+° 

^•"'-'-0.28 


1.17; 
5 42+2 '''' 
1.02; 


9.6 
2.1 
8.5 
5.7 


1.89 
0.37 
1.99 
1.07 


4.1 
3.2 
7.3 
7.1 


NGC 5408 


•^-0.88 


2.50;; 


4.1 


1.92 


4.08!2;« 


2.54;;; 


7.7 


1.30 


6.9 


IC5152 


2 35+2-2' 


1 63+"''* 


5.3 


1.17 


C QO+0.10 

J-^S_0.64 


C TQ+0.10 

3-i8_0.64 


12.2 


1.07 


3.9 



Notes: (1) the name of the galaxy; (2) the median of the quadratically-added amplitude of the non-circular motions (unconstrained case), averaged 
over the entire radial range, the errors indicate the lower and upper quartile; (3) same as (2), but averaged over the inner 1 kpc only; (4) the 
percentage the non-circular motions contribute to the maximum rotation velocity; (5) the median of the reswri absolute residual velocity field after 
the harmonic decomposition; (6), (7), (8), (9) the same as (2), (3), (4), (5), but for the constrained case; (10) maximum radius for the averaging of 
the quadratically-added amplitudes. 
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with po being the central density and the core radius. The rota- 
tion curve corresponding to this density profile can be expressed 
as 



Vrot('') = \j'^^Gpori 1 1 + — arctan 



(9) 



5.2. Fitting process 

As dwarf galaxies are believed to be dark matter dominated at 
all radii ( |de Blok & McGaugh|| 1997) 1, the velocity contribution 
from the baryons is often neglected when decomposing the ro- 
tation curves of these galaxies. But although dark matter is the 
dominant component in dwarf galaxies, the baryons, i.e., gas and 



stars, are still important, especially close to the centre. Including 
all components, the observed rotation velocity is then given by 



•"obs 



2 ,2,2 



halo' 



(10) 



We model the rotation curves in various ways, partly including 
the baryons. The different approaches are described in the fol- 
lowing. 

5.2.1. Minimum-disc case 

This is the simplest case ignoring all baryonic contributions. 
Thus, only the dark matter haloes are fitted to the observed rota- 
tion curve. 
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Fig. 9. Median of the absolute residual field from the harmonic 
decomposition (constrained case) vs. the one from the rotation 
curve analysis. The light grey triangles represent the values de- 
rived by CT08. They are complemented by the results of our 
analysis (black open triangles). The error bars indicate the lower 
and upper quartile. 

5.2.2. Minimum-disc -i- gas case 

In this case, the contribution of the gas is taken into account 
as well. Next to Hi, which dominates the gas component. He 
and metals are included by scaling the H i column density by a 
factor of 1.5. As the ratio between the molecular gas and Hi is 
much lower in dwarf galaxies than in luminous spirals (Taylor 
et al. 1998[ Leroy et al. 2005 1, we do not correct for this gas 
component. In order to create a surface density profile for the 
Hi, the GIPSY task elUnt is taken, using the Hi intensity maps 
(Figs.[2]to[7] lower right panels) and the parameters from the best 
tilted-ring model (see Table [3]l as an input. The output of ellint 
is given in terms of mean flux and has therefore to be converted 
to physical units. We then use the GIPSY task rotmod in order 
to determine the rotation of the gas under the assumption of an 
infinitesimally thin disc. 



5.2.3. Maximum-disc case 

Here, the contribution of the gas, the stars and the dark mat- 
ter halo are simultaneously fitted to the observed rotation curve. 
For all galaxies, the surface density can be described by a Sersic 
model ( |Sersic|1968"l l 



= /Jo 



1.086 



(11) 



with yUo being the central surface brightness, h the disc scale 
length, and n the shape parameter (where n - \ gives an ex- 
ponentialprofile). The parameters for each galaxy are listed 
in Table pi They are converted from mag/arcsec^ to units of 
Lo/pc^. Again, we use rotmod to determine the contribution of 
the stars under the assumption of an infinitesimally thin disc. The 
derivation of the contribution of the stellar disc is the most criti- 
cal part; first of all, the surface photometry depends on different 



factors which are difficult to estimate like, e.g., the extinction. 
Secondly, in order to derive the stellar rotation curve, the mass to 
light ratio M/L has to be known. Note that we make the assump- 
tion of a stellar M/L of 1 in the case of NGC 2366 and of 2 in 
the case of ESO215-G7009, which are in good agreement with 



adopting a Kroupa or a Kennicutt initial mass function ( Portinari 
|et al. '2004, and references therein). For the other sample galax- 
ies, calculations by assuming a Kroupa or Kennicutt initial mass 
function reveal stellar M/L ratios that result in a stellar rotation 
curve which lies above the observed one. Therefore, we choose 
smaller stellar M/L ratios, but always greater than 0.2, which is 

[2008) 1. In case of 



still a physical value (see also Spano et al. 



NGC 5408, we do not fit a stellar component as the M/L value 
has to be far below 0. 1 in order to fit the observed rotation curve. 

6. Results and discussion 

As described in, e.g., vE09a, the outer parts of the rotation curves 
show large uncertainties due to the sparsely filled tilted rings. 
Therefore, we decompose the rotation curves twice for each of 
the different cases (see Sect. |5.2| i and for each dark matter den- 
sity profile: first, using the entire rotation curve and secondly, 
using a truncated rotation curve restricted to the inner few kpcs 
only. The results for all sample galaxies are shown in Figs. [TO] 
to 15 From top to bottom the minimum-disc, the minimum- 
disc + gas and the maximum-disc cases are fitted. Solid trian- 
gles represent the observed rotation curve, long-dashed lines the 
dark matter halo component, dotted lines the gas component, 
and dashed-dotted lines the contribution of the stars. The result- 
ing model fit is indicated by a solid line. 

Tables|6]and|7]list the results of the best fits for both halo pro- 
files including the reduced values. We assume V200 to be the 
rotation speed, which we estimate from the flat parts of the ob- 
served rotation curves. Following Eq.|7) we then derive the con- 
centration parameter c, which is between 9 and 10 for all sample 
galaxies, i.e., consistent with CDM simulations. Both parame- 
ters are kept fixed so that r2oo is the only free parameter. In the 
case of the ISO halo both parameters, the core radius and the 
core density po, are allowed to vary. 

In the following subsections we present the results for each 
galaxy separately. 

6.1. NGC 2366 

A mass decomposition of NGC 2366 has already been done by 
Oh08. They derived the rotation curve from a so-called "bulk" 
velocity field, which only contains the circularly rotating com- 
ponents (Oh08). Our analysis is done on the Hermite velocity 
field that still contains non-circular motions. However, as shown 
in Sect. |4.1[ the non-circular motions in NGC 2366 (as well as in 
the other sample galaxies) do not significantly affect the rotation 
curves. 



Figure 10 shows the resulting model fits. A comparison of 
the reduced~jp values (see Tables |6] and [t]) reveals that the ISO 
halo gives much better results than the NFW halo independent of 
the radius of the rotation curve. In the outer parts, the observed 
rotation curve is affected by uncertainties caused by the sparsely 
filled tilted rings and by asymmetries in the H i velocity field (cf. 
vE09a). Therefore, it strongly declines from a radius of 5 kpc on 
so that we also model the curve within the inner 5 kpc only. This 
gives improved fits for both haloes: x^^^i is lower by a factor of 
four when using the NFW halo and even by a factor of 20 when 
using the ISO halo. Additionally, decreases by including gas 
and stars. 
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Table 5. Parameters for the exponential disc fitting. 



Parameter [Unit] 


NGC 2366 


ESO059-G001 


ESO215-G7009 


NGC 4861 


NGC 5408 


IC5152 




(1) 


(2) 


(3) 


(4) 


(5) 


(2) 


Band 


V 


H 


R 


R 


J 


H 


fio [mag/arcsec^] 


22.82 


20.23 


23.16 


22.95 


18.29 


17.65 


h [kpc] 


1.59 


0.81 


2.16 


0.88 


0.33 


0.36 


n 


1 


1.33 


1 


1 


1 


0.96 


MIL 


1 


0.8 


2 


0.3 


< 0.1 


0.2 



References: photometry by (l) |Hunter et al.| ( |200T) , (2) |Kirby et aTlpOOSl l, (3) |Warren et"aLlp006) . (4) |Gil de Paz & Madore|(2005| l, (5) |Noeske | 
[etaLlp003) . 

Table 6. The basic parameters of the mass decomposition - the NFW model. 









NFW full 






NFW small 






Parameter 


V200 


c 


'■200 


Arzoo X,^i 


V200 


c 


''200 


Arioo 




Unit 


[km s-'] 




[kpc] 


[kpc] 


[kms-'] 




[kpc] 


[kpc] 






(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 2366 


Min.-disc 


50 


9.4 


117.54 


7.26 47.47 


50 


9.4 


92.81 


4.52 


19.88 


Min.-disc + gas 


50 


9.4 


206.12 


16.44 51.76 


50 


9.4 


136.75 


5.06 


8.19 


Max. -disc 


50 


9.4 


386.29 


37.63 36.27 


50 


9.4 


229.04 


11.39 


7.81 


ESO059-G001 


Min.-disc 


69 


9.2 


89.22 


7.84 22.81 












Min.-disc + gas 


69 


9.2 


96.00 


8.04 19.67 












Max. -disc 


69 


9.2 


124.40 


21.27 63.70 












ESO215-G7009 


Min.-disc 


77 


9.0 


134.57 


5.13 19.07 


77 


9.0 


128.08 


3.63 


9.14 


Min.-disc + gas 


77 


9.0 


157.35 


7.85 30.18 


77 


9.0 


145.41 


5.00 


12.48 


Max. -disc 


77 


9.0 


309.86 


17.19 19.95 


77 


9.0 


283.81 


11.67 


9.13 


NGC 4861 


Min.-disc 


46 


9.5 


125.79 


10.55 17.42 


46 


9.5 


132.55 


18.92 


29.44 


Min.-disc + gas 


46 


9.5 


181.48 


13.27 10.93 


46 


9.5 


189.54 


24.43 


16.57 


Max. -disc 


46 


9.5 


219.91 


26.36 24.86 


46 


9.5 


267.28 


65.64 


37.31 


NGC 5408 


Min.-disc 


53 


9.7 


182.70 


24.22 47.97 


53 


9.7 


187.81 


34.65 


63.22 


Min.-disc + gas 


53 


9.7 


233.18 


30.95 38.76 


53 


9.7 


239.76 


46.08 


51.79 


IC5152 


Min.-disc 


44 


10.0 


34.26 


1.80 4.92 


44 


10.0 


31.03 


1.90 


4.51 


Min.-disc -I- gas 


44 


10.0 


43.47 


2.72 7.09 


44 


10.0 


37.09 


2.52 


5.12 


Max.-disc 


44 


10.0 


79.97 


5.57 6.93 


44 


10.0 


90.67 


13.45 


10.58 



Notes: (1) and(6): the rotation velocity at r2oo, estimated from the flat part of the rotation curve; (2) and (7): the concentration parameter derived 
following Eq.pj (3) and (8): the fitted virial radius; (4) and (9): uncertainties of ''200; (5) and (10): reduced values. 



The reduced H i rotation curve fitted with an ISO halo gives 
the lowest values xl^i- The resulting model fit agrees very 
well with the observed curve. Moreover, the values of the core 
densities and radii from the ISO models are comparable to the 
results of Oh08. 



6.2. ESO059-G001 

As the rotation velocities of this galaxy have small errors (also 
in the outer parts), we only perform a mass decomposition 



for the entire rotation curve. Figure 11 shows the resulting 
curves. Again, the ISO halo gives much better results than the 
NFW halo. Especially the models for the minimum-disc and 



minimum-disc + gas cases using the truncated rotation curve re- 
sult in very low;^'^^^ values (see Tables |6] and |7ji. 



6.3. ESO215-G7009 



ESO215-G7009 is modelled over the entire radial range and 
over the inner 8.1 kpc only (see Fig.[T2]|. In both cases, the NFW 



halo fits sometimes even give lower values than the ISO halo 
fits. This can probably be explained by the slightly scattered ro- 
tation velocities between 1 and 3 kpc. 
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NFW halo 



NFW halo < 5 kpc 



ISO halo 



SO halo < 5 kpc 
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_^lLJIIt*Ui'*'_ 


/ 


■ Min. Disc -1- Gas 








.1.1.1 


1,1,1,1,1.1.1,1, 
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Fig. 10. Mass decomposition of NGC 2366. The rotation curves for an NFW halo and a pseudo-isothermal halo are modelled for the 
full rotation curve, and for the inner 5 kpc only. From top to bottom the minimum-disc, the minimum-disc + gas, and the maximum- 
disc cases are fitted. Solid triangles represent the observed rotation curve, long-dashed lines the dark matter halo component, dotted 
lines the gas component, dashed-dotted lines the contribution of the stars, and solid lines the resulting model fit. 
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Fig. 11. Mass decomposition of ESO059-G001. The same as in Fig. 



10 



No reduced rotation curve is fitted. 
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Fig. 13. Mass decomposition of NGC4861. The same as in Fig. 



10 



A reduced rotation curve is fitted up to a radius of 5.5 kpc. 
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ISO full 








ISO small 








Parameter 




Po 


Apo 






Po 


Apo 




Ar, 


^red 


Unit 




[10-3Mopc-3] 


[lO-'Mopc-^] 


[kpc] 


[kpc] 


[10-3Mopc-3] 


[10-3Mopc-3] 


[kpc] 


[kpc] 








(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 2366 


Min.-disc 




56.66 


10.42 


0.94 


0.11 27.82 


32.95 


1.33 


1.58 


0.05 


1.79 


Min.-disc + 


gas 


53.12 


14.46 


0.77 


0.13 30.27 


27.54 


1.02 


1.44 


0.04 


0.88 


Max. -disc 




25.28 


9.82 


0.92 


0.22 24.99 


11.27 


0.72 


2.23 


0.14 


1.24 


ESO059-G001 


Min.-disc 




67.08 


3.18 


1.48 


0.06 0.84 












Min.-disc + 


gas 


62.98 


3.58 


1.49 


0.07 1.10 












Max. -disc 




25.70 


5.37 


3.03 


0.76 12.52 












ESO215-G7009 


Min.-disc 




77.74 


10.73 


1.27 


0.11 18.42 


71.67 


9.27 


1.35 


0.12 


15.67 


Min.-disc + 


gas 


98.86 


15.22 


1.01 


0.10 16.47 


89.34 


11.62 


1.10 


0.09 


11.77 


Max.-disc 




34.27 


8.49 


1.47 


0.24 19.18 


28.47 


6.10 


1.71 


0.26 


14.44 


NGC 4861 


Min.-disc 




13.98 


1.58 


2.22 


0.18 4.29 


11.06 


0.70 


3.04 


0.22 


1.00 


Min.-disc + 


gas 


9.67 


1.29 


2.45 


0.24 3.91 


8.75 


0.98 


2.90 


0.35 


1.72 


Max.-disc 




4.85 


1.13 


3.93 


0.81 10.40 


3.95 


1.48 


7.10 


7.19 


14.21 


NGC 5408 


Min.-disc 




9.22 


2.11 


3.36 


0.68 20.68 


6.69 


0.87 


6.37 


1.66 


6.85 


Min.-disc + 


gas 


6.98 


1.76 


3.70 


0.87 18.13 


4.94 


0.78 


8.64 


4.05 


7.24 


IC5152 


Min.-disc 




378.79 


40.76 


0.37 


0.02 1.94 


389.45 


40.23 


0.36 


0.02 


1.24 


Min.-disc + 


gas 


369.69 


48.16 


0.34 


0.03 2.23 


398.06 


49.18 


0.33 


0.03 


1.36 


Max.-disc 




42.11 


6.70 


1.18 


0.14 2.87 


47.24 


15.81 


1.08 


0.36 


4.44 



Notes: (1) and (6): fitted core density po of the pseudo-isothermal halo model; (2) and (7): uncertainties of po; (3) and (8): fitted core radius of 
the pseudo-isothermal halo model; (4) and (9): uncertainties of r^; (5) and (10): reduced values. 



6.4. NGC 4861 



6.6. IC5152 



Figure 13 shows the results for NGC 4861. The reduced val- 
ues of the maximum-disc case reveal that the stars are not well 
represented by our estimate as the x^^^i values increase, inde- 
pendent of the model, after adding the stellar component (see 
Tables [6] and [Tji. Similar to the other galaxies in our sample the 
ISO matches the observed rotation curve better than the NFW 
halo. 



6.5. NGC 5408 

The H I rotation curve of NGC 5408 rises very slowly in the in- 
ner 2 kpc. Using the photometric data from Noeske et al. ( 2003| l 
(see Table |5]), it is impossible to derive a stellar rotation curve 
that lies below the total rotation curve, even when assuming a 
stellar M/L ratio of 0. 1, which is the lowest value before becom- 
ing unphysical ( [Spano et al.|2008| . We either underestimate the 
rotation in the central part of the galaxy due to the low spatial 
resolution or the photometric data are not accurate enough. We 
therefore decide against modelling the stellar component. 

shows the resulting fits. Again, the;^'^^^^ values are 



Figure 
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very high for the NFW halo and much smaller for the ISO halo. 
For this galaxy, the reduced rotation curve did not improve the 
quality of the models, which is probably due to the fact that 
the repeated rise of the rotation curve beyond a radius of 7 kpc 
matches the first rise of the curve. 



The results for IC5152 are presented in Fig. 15 We fit NFW 
and ISO haloes to the entire radial range and to the inner 2. 1 kpc 
only. Although both dark matter halo models provide adequate 
fits (the x^^^i values for the NFW model are the lowest of our 
sample), the ISO model still gives better values than the NFW 
halo. 



6.7. Discussion 

In contrast to recent work by dB08 and Oh08, we fitted the NFW 
profiles with fixed values of V20() and c (see the beginning of 
this section). With this approach we assumed reasonable values 
for V200 and c that follow cosmological predictions, and also got 
reasonable results for r2oo. Nevertheless, we could show that in 
general, the observed rotation curves can better be fitted by an 
ISO halo than by an NFW halo, although at least in the case of 
ESO215-G7009 the NFW fits are equally good or even better 
than the ISO fits. The values of the core radius rc and the core 
density po are plausible and in good agreement with results for 
LSB galaxies (de Blok et al..2001a, Kiizio de Naray et al.|2008) l 
and for dwarf galaxies ( Spario et al.|2008 1. They also follow the 
recent finding of a constant central surface density //qd - '"cPo 
(Donato et al. 2009 and references therein). 

As mentioned in the introduction, CDM simulations predict 
density profiles with a steep inner slope of a = -1 ( |Navarro| 
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Fig. 14. Mass decomposition of NGC 5408. The stellar contribution cannot be fitted. A reduced rotation curve is fitted up to a radius 
of 6.3 kpc. Otherwise the same as in Fig. 10 
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Fig. 15. Mass decomposition of IC 5 152. The same as in Fig. 10 A reduced rotation curve is fitted up to a radius of 2. 1 kpc. 



let al. 



T9961 ) or even a = -1.5 ( |Moore et al.|T998l l. |de Blok et al.| 



( 200 a I showed for their sample of LSB galaxies that the ob- 
served distribution of a peaks at -0.2 + 0.2. In order to compare 
these results with ours, we first measure the slopes in the inner 
kpc of our data by plotting the logarithm of the density vs. the 



logarithm of the radius (see Fig. 16 1. For all galaxies except for 
1C5152, a lies between -0.43 and 0.03 (see also Table[8)l, which 
is within the errors of the a determined by |de Blok et al.| ( )2001a^ . 



ESO215-G7009, which partly showed higher values for the 
ISO halo model than for the NFW halo model (see Sect. |6.3| l, 
has an a of -0.35. This means that its innermost slope can better 
be described by an ISO halo. 

In Fig. [16] we also fitted the ISO and NFW profiles of the 
minimum-disc case (dotted and long-dashed lines respectively). 
They do not always agree with the observed density profiles, 
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Galaxy 


a 


Aa 


'■in 








[kpc] 




(1) 


(2) 


(3) 


NGC 2366 


0.03 


0.08 


0.10 


ESO059-G001 


-0.29 




0.52 


ESO215-G7009 


-0.35 


0.23 


0.28 


NGC 4861 


-0.01 




0.50 


NGC 5408 


-0.43 




0.35 


IC5152 


-0.82 


0.07 


0.23 



Notes: (1) the slope of the inner 1 kpc of the density profile; (2) uncer- 
tainties; (3) the radius of the innermost point. 



which is to be expected as x^^.^ is often quite high indicating a 
bad fit. 



We now take the plot by |de Blok et al. ( 2001a i (see also 
Oh08), which shows a vs. the logarithm of the innermost radius 
(see Fig. 17 1. Data from previous studies are plotted with grey 
open circles ( de Blok et al.|200Tal l, grey open squares (de Blok 



& Bosma_200 2|l, grey open triangles ( |Swaters et al. 2003| ), and 
black solid triangles (Oh08). Our own results are plotted with 
black solid squares and are in good agreement with all previous 
works. Only IC 5152, which has an inner slope of a = -0.82, is 
close to CDM predictions. 

We now discuss possible systematic errors in the data 
(as mentioned in the introduction): Swat ers et al.| ( |2003| ) and 
Spekkens et al.| (j2005 ) used rotation curves of dwarf and low 
surface brightness galaxies obtained from optical long-slit spec- 
tra in order to measure the slopes of the mass distribution. 
They compared their results to model spectra by taking into ac- 
count possible sources of uncertainties like slit misalignment, 
slit width, and seeing. Both studies came to the result that the 
galaxies do not require haloes with steep cusps, but that haloes 
with Q- = 1 cannot be ruled out. However, they clearly showed 
that a steeper slope (a > 1) is not consistent with the observed 
profiles. As we use 2d velocity fields, our data are not affected 
by the above mentioned errors. We can also rule out geometric 
effects caused by low velocity gas picked up along the line of 
sight at higher inclinations as our data have a very high spectral 
resolution. 



As already mentioned in Sect. 4. 1 CT08 could show that in- 
accurate centre positions do not change the inner slope of a rota- 
tion curve. Non-circular motions are detected, but small enough 
to not significantly change the rotation velocities. Nevertheless, 
there are problems that affect our results as, e.g., the low spatial 
resolution of especially the ATCA data. 

Summarised it can be said that for the majority of galax- 
ies studied here, the ISO halo reproduces the observed rotation 
curves much better than the NFW halo, which is in good agree- 
ment with previous studies (e.g., |Kuzio de Naray et al. 2008 



|Spano et al. 2008 ). The slopes of the observed rotation curves 
are inconsistent with CDM predictions, but agree with the results 
from de Blok et al. ( 2001a) l. One explanation for the cusp-core 
discrepancy might be that most simulations neglect the baryons. 
And indeed, first approaches to include baryons into CDM sim- 
ulations show that their contribution changes the inner density 
profile from a cusp to a core ( Romano-Diaz et al.|2008 1. 




-1 
log r,^ (kpc) 



Fig. 17. The inner slope of the dark matter density profiles plot- 
ted against the radius of the innermost point. Presented in grey 
are the results from publications by de Blok et al. ( 2001a i 
(open circles), de Blok & Bosmaj(2002|l (open squares), and by 
|S waters et al.| ( |2003[ ) (open triangles). Shown in black are the 
results from Oh08 (solid triangles). Our results are overplotted 
with black squares. We fit the isothermal and NFW profiles fol- 
lowing [deBEketig (|200g . 



7. Summary 

We used VLA and ATCA Hi synthesis data of a sample of 
six nearby irregular dwarf galaxies in order to decompose the 
observed rotation curves into the contributions from the stars, 
the gas, and the dark matter halo. In order to rule out system- 
atic effects, we first performed a harmonic decomposition of the 
Hermite velocity fields. The results show that the quadratically- 
added amplitude Ar(r) is generally below 6kms"', in the cen- 
tral kpc even below 3kms"^. Over the entire radial range non- 
circular motions contribute less than 25%, sometimes even 
less than 10% to the local rotation velocity (constrained case). 
According to [Hayashi et al. ( |2004 1, the triaxial dark matter 



haloes cause non-circular motions as large as 50% of the lo- 
cal rotation velocity within the central kpc. Our much smaller 
amplitudes demonstrate that the observed rotation curves of the 
galaxies studied here are not severely influenced by non-circular 
motions. 

With this knowledge we decomposed the observed rotation 
curves into their components. We fitted a cuspy NFW halo 
and a cored pseudo-isothermal halo to the rotation curves and 
modelled the minimum-disc, the minimum-disc + gas, and the 
maximum-disc cases. The resulting model fits show that in all 
sample galaxies (except for ESO215-G7009), the observed ro- 
tation curves are better represented by the ISO halo than by the 
NFW halo. As the outer parts of the rotation curves are affected 
by uncertainties caused by the sparsely filled tilted rings, we also 
performed a mass decomposition by fitting the inner few kpcs 
only. For most galaxies, this significantly improved the qual- 
ity of both the NFW and the ISO halo models. The reduced 

could in most cases be decreased for both halo models by 
including the gas. In some cases, a further improvement could 
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Fig. 16. The density profiles of all sample galaxies derived from the observed rotation curves (open grey triangles). Their inner 
slopes a are measured by applying a least square fit to all data points within the innermost kpc (bold black lines). The fitted values 
of a and the uncertainties are placed into the upper right corner of each panel. Note that the rotation curves of ESO059-G001, 
NGC 4861, and NGC 5408 only contain two points in the inner 1 kpc. Therefore, no uncertainties can be given. The long-dashed 
and dotted lines show the NFW and the ISO profiles respectively using the parameters of the minimum-disc case. 



be achieved when also including the stars. This implies that the 
baryons should not be neglected, at least not in the inner part of 



dwarf galaxies (see also S waters et al.|2009) 
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